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Abstract 

We present theoretically the electronic structure of antiferromagnetic (AFM) FeSe monolayer on 
TiC>2 terminated SrTiO3(001) surface. It is revealed that the striking disappearance of the Fermi 
surface around the Brillouin zone (BZ) center can be well explained by the antiferromatnetic (AFM) 
phase. We show that the system has a considerable charge transfer from SrTiO3(001) substrate to 
FeSe monolayer, and so has a self-constructed electric field. The FeSe monolayer band structure 
near the BZ center is sensitive to charge doping, and the spin-resolved energy bands at BZ corner 
are distorted to be flattened by the perpendicular electric field. We propose a tight-binding model 
Hamiltonian to take these key factors into account. We also show that this composite structure is an 
ideal electron-hole bilayer system, with electrons and holes respectively formed in FeSe monolayer 
and TiC>2 surface layer. 
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Due to its simple crystal structure, prominent antiferromagnetic magnetic ordering, and 
significant pressure effect on superconductivity, during the last few years FeSe binary su- 
perconductor has been drawing extensive attention from fundamental mechanism of Fe- 
based superconductivity to promising application in such areas as superconducting wires 
and thin films [l]. Recently, monolayer FeSe with atomic flatness was successfully grown on 
SrTiO3(001) surface 2], and strikingly, the high-temperature superconductivity signature 
from the monolayer was revealed by high energy resolution scanning tunneling spectroscopy 
(STS). It is even hopeful that with further improving the sample quality, the transition 
temperature is probably as high as 77 K, the liquid nitrogen boiling point. The subsequent 
angle-resolved photoemission spectroscopy (ARPES) experiment [3| shows that the Fermi 
surface of this FeSe monolayer only consists of electron pockets near the Brillouin zone (BZ) 
corner, with no hole pockets around the zone center observed. Theoretical efforts are being 
paid to explain these intriguing STS and ARPES reports. The soft phonon of SrTiO 3 (001) 
was found to strengthen the Cooper pairing {4]. £4 symmetry was also noticed [5|, wherein 
it was suggested that cuprates and iron-based superconductors share an identical high-T c 
superconducting mechanism. The density functional theory (DFT) total-energy calculations 
gave that collinear antiferromagnetic order (CAFM) FeSe monolayer on Ti02 terminated 
SrTiO3(001) surface is the most stable structure, and there is neither hybridization nor 
charge transfer between FeSe monolayer and SrTiC^OOl) surface [6]. The calculated Fermi 
surface for this CAFM structure exists around the M point in the fold BZ, and there is a 
Dirac-cone-like bands of FeSe near Fermi energy, which unfortunately are not observed by the 
ARPES measurement [3| that the Fermi surface disappears around the zone center. A more 

Q 

recent DFT study [7J on a single FeSe monolayer (without substrate) showed that an AFM 
instead of CAFM phase can have the Fermi surface resembling the ARPES data. Basically, 
however, the interaction between AFM-ordered FeSe monolayer and SrTiO3(001) surface re- 
mains largely unclear in literature, despite its urgent importance on understanding the key 
role played by the interface in controlling the Fermi surface topology and high-temperature 
superconductivity of FeSe monolayer. 

In this work, we carefully investigate with DFT based ab initio calculations the elec- 
tronic structure of AFM FeSe monolayer on Ti02 terminated SrTiO3(001) surface. We 
show that the calculated electronic structure of the system is close to the ARPES experi- 
mental results. The charge transfer and electric field are induced by the interaction between 
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SrTiO3(001) substrate and AFM FeSe monolayer, and turn back out to play an important 
role in modulating the band structure and Fermi-surface topology of the FeSe monolayer. 
One extraordinary consequence is that the spin-resolved electron pockets surrounding the BZ 
corner are distorted and become flattened. We propose a tight-binding model Hamiltonian 
that can describe these key factors. 

We use 6-layer SrTiOs (001) slabs to mimic the surface of SrTiOs substrate, as shown in 
Fig. Ufa). The SrTiOs(OOl) surface is either SrO or Ti02 terminated with the nearly same 
surface energy js-10|. According to the experiment we use the Ti02 terminated surface 
in our atomic models. FeSe monolayer is adsorbed on the lxl two-dimensional unit cell of 
TiCVterminated SrTiOs (001) surface. To test the influence of SrTiOs antiferrodistortion on 
surface electronic structures, we also calculated a y/2 x \/2 supercell, in which the lower four 



layers of Ti0 2 and SrO are fixed to t 
is set to the experimental value 2.1° 



leir bulk positions and the octahedral rotation angle 



11] . In our models, the vacuum space is larger than 



10 A. The results are also checked for vacuum space larger than 20 A. The bulk calculation 
shows that the lattice parameter of SrTiOs is 3.901 A, which agrees well with experiments 
and previous theoretical results. The adsorption structure relaxation is performed with fixed 
lattice parameters. 

The total energy and electronic structure calculations are performed by using projector 



augmented wave (PAW) method 



14| . The exchange correlation potential is describe d by 



the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) type 15]. 
The plane wave cut-off energy is chosen to be 400 eV, which is converged in our test. 9x9x1 
and 6x6x1 Monkhorst-Pack fc-points are used in the reciprocal space integration for lxl 
and \/2 x \/2 supercells. The system is relaxed until the force on each atom is smaller than 
0.01 eV/A. To include the strong correlation effect for d electrons of Fe atoms, we adopt 
Hubbard-[7 correction. We use C/=0.5 eV for Fe-3d electrons in our main calculations, and 
the influence of different values of Hubbard-[7 on electronic structures are also discussed 
in the following text. All the DFT calculations are performed by using Vienna Ab-initio 
Simulation Package (VASP) jl2]. The relaxed adsorption structure of FeSe monolayer on 
Ti02 terminated SrTiOs (001) surface is shown in Fig. [TJ The lower layer Se atoms are on 
the top of Ti atoms, while the Fe atoms locate on the top of O atoms. The vertical distance 
between Ti and Se atoms is 3.13 A, which is slightly larger than that of the CAFM phase by 
0.07 A. The vertical distance between Fe and O atoms is 4.43 A. The O atoms in the Ti02 
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FIG. 1: (Color online) The side view (a) and top view (b) for the atomic structure of monolayer 
FeSe absorbed on Ti02 terminated SrTiC^OOl) surface. The magnetic ordering is also represented 
in top view (b) as red arrows. 

top-layer arise slightly, they are 0.1 A higher than Ti atoms. The magnetic moment on each 
atom is calculated by integrating the spin density in Wigner-Seitz radius of 1.30 A, 1.16 A, 
2.14 A, 1.22 A, and 0.82 A for Fe, Se, Sr, Ti, and O atoms respectively. The calculated 
magnetic moment is 2.4 /ig on each Fe atom and is negligible on the other atoms. 

The energy bands for spin- up channel are shown in Fig. [2(a). The spin-down energy 
bands are the same as those for spin-up channel except for an exchange of X and X' points 
in BZ. There are two bands [indicated as a and 7 in Fig. 12(a)] crossing the Fermi level 
Ep near M point and foming an ideal pair of electron-hole bands with perfect Fermi-surface 
nesting, while around Y point no metallic bands occur. Thus, the calculated Fermi surface 
only exists at the BZ corner, which agrees with the ARPES result in experiment 3j . In our 
calculations, the antiferrodistortion of SrTiOs substrate does not have considerable influence 
to the electronic structure of adsorbed AFM FeSe monolayer. 

The disappearance of Fermi surface at BZ center is the result of a full electron filling of 
band (3. We find that the shape of band j3 can be altered dramatically by changing the value 
of Hubbard-?/ correction. Actually, band shall cross the Fermi level at the BZ center in 
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FIG. 2: (Color online) The spin-up energy bands of monolayer FeSe adsorbed on TiO"2 terminated 
SrTiO^OOl) surface (a). The red and blue lines are contributed by 3d orbitals on Fe atoms and 2p 
orbitals on O atoms, respectively. The Fermi level is set to be zero. Panel (b) shows the influence 
of Hubbard- U on the energy difference Ap. 

the absence of Hubbard-?/ correction. When the Hubbard-?/ correction is applied, band (3 
lowers its energy and a gap of 41 meV [indicated by green arrows in Fig. [2(a)] is obtained for 
£7=0.5 eV, which agrees with the experimental result 3j. The detailed relationship between 
the minimum energy difference of band and Ep around the BZ center (denoted as Ap) and 
the Hubbard-t/ correction is shown in Fig. 12(b). One can see that Ap and U are mainly in 
a linear relation. We have checked that this phenomenon does not exist in the paramagnetic 
FeSe/SrTiC>3 system, wherein there are energy bands crossing the Fermi level around the 
BZ center even for 5.0 eV Hubbard-?/ correction. 

Our PDOS (partial density of states) analysis shows that bands a and /3 are contributed 
by Fe 3d orbitals, while the hole-type band 7 comes from the O 2p orbitals. We analyze the 
wavefunctions of band /3 at T point and bands a, (3 and 7 at M point. The wavefunction 
of band /3 at T point shows a shape of a dumbbell in z-direction with a torus in xy-plane, 
which is the typical d z i feature. Whereas, the wavefunctions of bands a and (3 at M point 
show vertical four-leaf shapes, which are the d xz and d yz features. The wavefunctions for 
band 7 are dumbbell shaped p x and p y orbitals that are mainly distributed on the surface 
O atoms. The crossing between bands a and 7 indicates that charge transfer occurs from 
surface O atoms of SrTi03 substrate to Fe atoms of FeSe monolayer. Our Bader analysis 
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identifies a charge transfer of —0.17 e per unit cell. 

From above results, we know that the AFM FeSe monolayer is charge doped by the SrTiC>3 
substrate. In order to analyze this effect, we study a free-standing AFM FeSe monolayer by 
doping charge from —0.5 e to 0.5 e. The results are shown in Fig. [3J We see that the shape 
of band /3 changes dramatically when the system is positively charged. Band at V point 
crosses the Fermi level, and becomes a hole pocket. When the FeSe monolayer is negatively 
charged, the shape of band keeps static, while band a at T point lowers its energy and 
finally crosses the Fermi energy, showing an electron pocket. Both positive and negative 
charge doping increase Ep relative to the energy of band a and at M point. These results 
agree well with the previous study on FeSe monolayer without Hubbard-f/ correction 
However, as we have determined above, the realistic charge transfer in FeSe/SrTiC>3 is only 
— 0.17e. By comparison with Fig. |3l we find that this value is not large enough for the 
existence of electron pocket at V point. Thus, the energy bands of AFM FeSe monolayer are 
altered limitedly by pure charge doping from SrTiC>3 substrate. 

Besides the charge doping, another prominnet feature in Fig. |2]is that bands a and f3 are 
nonsymmetric along M-X and M-X' directions, and a small splitting appears at M point. 
Whereas, for the free-standing AFM FeSe monolayer, bands a and /3 are symmetric along 
M-X and M-X' directions, and they are degenerate at M point for both neutral and charged 
systems as shown in Fig. |3j Considering the net charge and dipole effect of SrTi0 3 substrate, 
a vertical electric field may exist in FeSe/SrTi03. Then we further study the electric filed 
effect. In our further DFT calculations, a free-standing monolayer FeSe is particularly 
designed to be exposed in a perpendicular electric filed, and the atomic positions are relaxed 
for each electric field strength. The calculation results are shown in Fig. HI Clearly, when 
the electric field strength is zero, the energy bands are symmetric along M-X and M-X' 
directions, and the bands of spin-up and spin-down components are the same. When the 
electric field strength is nonzero, however, we are astonished to see that the energy bands 
along M-X and M-X' directions become nonsymmetric, with an obvious spin splitting at 
M point. The splitting amplitude increases with increasing the electric field strength. So 
we verify that the energy bands near the Fermi energy and M point are distorted by the 
perpendicular electric field, and the energy of band a at X and X' points are different. 
Besides that, the features of spin-up and spin-down bands are just opposite, namely, the 
energy bands in M-X' (M-X) direction for spin-up electrons are just the same as the energy 
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FIG. 3: (Color online) The electronic energy bands of free-standing monolayer FeSe with different 
levels of charge doping. The value of charge doping is changed from —0.5 e to 0.5 e per unit cell. 
The Fermi energy is set to zero. 

bands in M-X (M-X') direction for spin-down electrons. Comparing the spin-up bands in 
Fig. H] for free-standing FeSe monolayer and those for FeSe monolayer on SiTiOs surface 
(Fig. |2]), we see that the distortion of the energy bands are quite similar. Thus, the effect 
of self-established electric field in FeSe/SiTi03 is confirmed. 

To explore the detailed mechanism of energy band distortion, we construct a tight-binding 
model Hamiltonian to describe the system. As we have discussed above, the corresponding 
wavefunctions of bands a and near M point are composed of d xz and d yz orbitals of Fe 
atom, and the wavefunction of energy band /3 near T point mainly comes from d z 2 orbital 
of Fe atom. From Fig. H] we know that the energy band distortion mainly occurs near M 
point. Thus, we will use d xz and d yz orbitals of Fe atoms in our tight-binding model, and 
neglect the d z 2 orbital for simplicity. However, the vertical electric field could not affect a 
model which only consists of d xz and d yz orbitals of Fe atoms, since all the centers of these 
orbitals are in a horizontal plane. Thus some out-plane orbitals should be included in our 
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FIG. 4: (Color online) The electronic energy bands of free-standing monolayer FeSe with different 
vertical electric field strength. The value of electric field strength is changed from 0.0 eV/A to 1.0 
eV/A. The energy bands near Fermi energy are distorted and show different behaviors for spin- up 
and spin-down channels. The Fermi energy is set to zero. 

model. The most natural choice is to use the orbitals of top and bottom Se atoms. We 
assume the Se orbitals in our model are symmetric in xy-plane, which may be Se 4p z orbital 
or its combination with Se 4s orbital. The following calculation shows that this assumption 
well describes the energy bands along M-X and M-X' directions. 

In our mean-field tight-binding model of AFM FeSe monolayer, we describe the spin-up 
and spin-down electrons separately. For spin-up electrons, the tight-binding model contains 
d xz and d yz orbitals of the spin-up Fe atom, besides one symmetric orbital for each Se atom. 
Similarly, for spin-down electrons, the tight-binding model contains d xz and d yz orbitals of 
the spin-down Fe atom, and one symmetric orbital for each Se atom. The hopping of spin- 
up electrons is schematically shown in Fig. [5j For the spin-down electrons, whereas, the 
hopping schematic diagram can be obtained by switching the positions of top and bottom 
Se atoms. 
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FIG. 5: (Color online) The hopping schematic diagram of spin-up tight-binding model in a unit 
cell. The small red circles, black large circles and gray large circles show the positions of Fe atoms, 
top Se atoms and bottom Se atoms separately. The spin-down Fe atom is located at the cell center 
and is not plotted in this spin-up tight-binding model. The spin-down tight-binding model can be 
obtained by exchanging the positions of black and gray circles. 

Considering the symmetry of each orbital, the tight-binding model Hamiltonian for spin- 
up electrons is written as H ( up ) = H^^te + , where 



^onSite = Yl \ € <i( C x~z,r£xz,r + C+ z ^ yz f) + (e s + dE) C + p ^ top f + (e a - dE) C + t -C bot , 



r+x) "F ^2\C XZ fC X z,f+y ~F C,,^ ^C, 



yz,r^yz,r+xj 



(1) 

(2) 



+h(ctz,'i iC bot,r+x/2 + Cy Z<i fitop,r+y/2) 



+ h.c. 



Here, c xz - and c+ - are annihilation operators of Fe d xz and d yz orbitals, while cf on p and 



yz,r 



top,r 



c totr are those of top and bottom Se symmetric orbitals, respectively. For simplicity the 
spin indices are omitted. The parameter d is the distance between top (or bottom) Se atom 
and Fe plane, and E are the electric field strength. 

The tight-binding model Hamiltonian for spin-down electrons is similar to that of spin-up 
electrons, except for that the positions of top and bottom Se atoms are exchanged. Thus, 



the Hamiltonian for spin-down electrons is written as ff( down ) = H^^tc + H' 



j-(down) 



where 



H 



(down) 
on-site 



(3) 



C yz,r C V z y+ £ ) 



(4) 



+t3{Cx zf £top,r+x/2 + C y Z ,f£bot,r+y/2) 



h.c. 



By performing the Fourier transformation, the model Hamiltonians for spin-up and spin- 
down electrons in momentum space are given by 



#< up) (k) 



e^(k) -it 3 sm(k y /2) 

e d Jk) 

it 3 sm(k y /2) e top 

it 3 sm(k x /2) 





-it 3 sin(/c x /2) 


Cbot 



(5) 



# {down) (k) 



(6) 



e dxz (k) -it 3 sm(k x /2) 

e dyz (k) -i* 3 sin(V 2 ) 

i^sin^^) e top 

it 3 sm(k y /2) e bot 

where e dxz (k)=e d +t 1 cos(k y )+t 2 cos(k x ), e dyz (k)=e d +t 1 cos(k x )+t 2 cos(k y ), e top =e s +dE, and 
ebot=e s — dE. 

After diagonalization of Eqs. (5) and (6), the FeSe monolayer's energy bands near the 
Fermi surface are plotted in Fig. El in which we have used the parameters of 6^=0.36 eV, 
e s =-4.57 eV, d=1.37 A, ^=-0.01 eV, t 2 =0.52 eV and t 3 =0.71 eV. The results show the 
Fermi surface distortion clearly. The energy bands are symmetric along X-M-X', and are 
the same for spin-up and spin-down electrons. In the presence of a nonzero vertical electric 
field, an energy splitting appears at M point. The amplitude of this splitting increases with 
increasing the electric field strength. The energy bands become asymmetric along X-M-X'. 
These features are quantitatively consistent with our DFT results [see Fig. (J4j)]. Thus, our 
tight-binding model further explains the origin of energy bands distortion in FeSe/SrTiCV 
Here, it should be noticed that in order to describe the electronic structures for the entire 
Brillouin zone, the tight-binding model should contains Fe d z 2 orbital besides the orbitals 
discussed above, since the energy band f3 in Fig. |2] at V point is mainly contributed by d z a 
orbital. 

To summarize, by performing the ab initio calculations, we have systematically inves- 
tigated the electronic structure of AFM FeSe monolayer on Ti0 2 terminated SrTiC^OOl) 
surface. We have shown that the interaction between SrTiO 3 (001) substrate and AFM FeSe 
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FIG. 6: (Color online) Calculated tight-binding energy bands of monolayer FeSe for spin-up (a-c) 
and spin-down (d-f) channels. The energy bands distortion is shown near M point. The Fermi 
energy is set to zero. 

monolayer induces charge transfer and a self-constructed electric field, which fundamentally 
modify the band structure and Fermi surface topology of the FeSe monolayer. In particular, 
the energy bands near the Fermi level are dramatically distorted by the electric field and 
show an obvious spin-resolved splitting around M point in BZ. As a result, the spin-resolved 
electron pockets evolved from band a become anisotropic and flat. We have proposed a 
tight-binding model Hamiltonian to reproduce these features. Also, we have shown that the 
present monolayer FeSe/SrTiO 3 (001) composite structure is an ideal electron-hole bilayer 
system, with electrons and holes respectively formed in FeSe monolayer and TiC>2. These 
findings are to be verified by the future spin-polarized ARPES experiments and are expected 
to shed light on the high-temperature superconductivity at the oxide interface. 
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